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ABSTRACT 

A flexible framework for the development of block structured volume grids for hypersonic Navier-Stokes 
flow simulations was developed for the analysis of the Shuttle Orbiter Columbia. The development of the 
flexible framework, resulted in an ability to quickly generate meshes to directly correlate solutions 
contributed by participating groups on a common surface mesh, providing confidence for the extension of 
the envelope of solutions and damage scenarios. The framework draws on the experience of NASA Langley 
and NASA Ames Research Centers in structured grid generation, and consists of a grid generation process 
that is implemented through a division of responsibilities. The nominal division of labor consisted of NASA 
Johnson Space Center coordinating the damage scenarios to be analyzed by the Aerothermodynamics 
Columbia Accident Investigation (CAI) team, Ames developing the surface grids that described the 
computational volume about the orbiter, and Langley improving grid quality of Ames generated data and 
constructing the final volume grids. Distributing the work among the participants in the 
Aerothermodynamics CAI team resulted in significantly less time required to construct complete meshes 
than possible by any individual participant. The approach demonstrated that the One-NASA grid 
generation team could sustain the demand for new meshes to explore new damage scenarios within an 
aggressive timeline. 


INTRODUCTION 

On Sunday February 1st, 2003, Space Shuttle Orbiter Columbia accident occurred during Earth 
re-entry while performing the Shuttle Transportation System mission 107 (STS-107). In response, the 
Columbia Accident Investigation (CAI) team was formed by NASA. Within the CAI team, the 
Aerothermodynamic CAI (ACAI) team investigated the role of hypersonic surface heat transfer on the 
demise if the Shuttle Orbiter. The flow simulation software that was used to assess the viscous hypersonic 
flow fields require block structured grids. These grids are difficult to generate because of the geometric 
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complexity of the flow field, and the vehicle shape, and the different mesh quality and mesh topology 
demands of the various flow simulation codes. 

This paper lays out the process that evolved into a flexible framework for the development of a 
common structured grid model of the Space Shuttle Orbiter, as needed by the ACAI. The cooperation 
among ACAI team members led to efficient grid generation through the separation of steps throughout the 
process. The result of this framework was the generation of a family of five volume grids in a short time 
frame, enabling the computational community to respond to the need for simulating viscous flow fields. 


A GRID GENERATION FRAMEWORK 

The flexible framework used by the ACAI for constructing grids consisted of a grid generation process and 
a division of responsibilities. Flexibility was enabled by permitting ACAI team participants to contribute 
at any level of the process. During the investigation, several participants performed grid quality analyses to 
ensure the viability of generated grids for their respective CFD tools. In order to develop a grid generation 
process for constructing structured grids about the Space Shuttle Orbiter, several issues required 
addressing: 

1. Flow solvers in use; 

2. topologies and grid resolution; and 

3. grid quality. 

Each issue is detailed in the following sections. Then, the process that was developed to address each issue 
is explained along with the rationale. Due to the fast pace of the investigation, the process chosen is very 
specific to the needs of the ACAI team, while not necessarily constituting the best practices of grid 
generation. Rather, the process is a compromise among all the issues in the above list, given the demands 
to support the availability of distributed grid generation resources, the need to have grid independent 
solution validation between flow solvers, and a very aggressive time line. 

FLOW SOLVERS 


Prior to constructing any discrete representation of a flow field domain, the requirements for the flow 
solvers in question needed to be addressed. For the case of the ACAI team, the flow solvers needed to 
accurately predict low density, viscous flow field phenomenon. The intended flow solvers were: the Langley 
Aerothermodynamic Upwind Relaxation Algorithm[l] (LAURA), the General Aerodynamic Simulation 
Program[2] (GASP), the Data Parallel Lower Upper Relaxation 3D (DPLUR3D) code, the Unified Solution 
Algorithm[3] (USA), and the Sandia Advanced Code for Compressible Aerothermodynamics Research and 
Analysis[4] (SACCARA). Each code has the ability to simulate the flow field phenomenon in question. 
However, each code has unique demands on the mesh. 

For example, the DPLUR3D requires that the block structured mesh have both face-to-face matching 
with no sub-facing and point-to-point matching across each block interface. The LAURA code, meanwhile 
can use a point-to-point, but can also have multiple cells matching a single cell on either side of the 
interface. Further, the LAURA code makes use of a grid adaptation scheme which requires that grid lines 
emanating from the OML surface traverse through the flow field and map directly to the outer boundary of 
the flow domain. Most flow solvers typically yield improved solution accuracy when orthogonality of the 
mesh to the OML is enforced. For good heat transfer resolution the spacing at and near the OML surface 
in the normal direction must be strictly enforced such that the first point off the surface is placed at a 
distance that is proportional to the local speed of sound. Near wall spacing and the number of points and 
distribution of the points within the boundary layer must meet flow solver accuracy requirements. Because 
some flow solver codes have difficulties with singularities, the ACAI team decided not to utilize pole 
boundaries in the common volume grids. Finally, grid lines skewness should be limited within the volume 
grid. 

Hence, the requirements for the grids to be produced are a direct result of the chosen flow simulation 
software. In particular, the volume grids should be orthogonal to the wall, have minimized skewness where 
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possible, without singularities, and provide a point-to-point mapping from the wall to the outer boundary 
and block faces. By maintaining these requirements on grid construction, the full suite of flow solvers 
mentioned should be able to provide accurate aerothermodynamic predictions of the hypersonic flow field 
phenomenon. 


TOPOLOGY AND GRID RESOLUTION 


Structured grids are characterized by three mutually orthogonal computational coordinates, forming a 
rectilinear lattice of points in an ordered three-dimensional (3D) array. When the mesh is combined with a 
mapping that dictates how the points are conformed to the surface of, and volume surrounding, the vehicle, 
a topology is created. The mapping is the link between the physical domain where the vehicle is described 
with physical coordinates such as (x,y,z) in a Cartesian domain, to a computational domain where the 
coordinates are typically (i,j,k), and typically form cube cells. Because of the requirements imposed by the 
various flow solvers and the need to develop a common mesh that could accommodate all the methods, the 
choice of mesh topology, structure, and quality were limited. As such, a multiple block grid with a global 
O-grid topology over the vehicle was chosen, coupled with a C-grid on the tail and wing, and a region along 
the wing leading edge that could accommodate an embedded O-grid, as shown in figure 1. 



Figure 1: Shuttle Orbiter structured grid topology. 


Shaded in yellow is the region of damage scenarios where the embedded O-grid could be used to isolate 
various types of geometry changes. The topological structures in use serve to limit the grid resolution 
where it is not warranted, while offering enhanced resolution were it is needed for flow solver accuracy. 

There are three basic possible structured grid topologies, as shown in figure 2. 

Each topology is so named because of its structure; the O-grid looks like an O, and similarly for the C- 
and H-grids. Each topology has advantages and disadvantages. For the shuttle geometry, a full body 
O-grid is placed around the fuselage with 235 points in the streamwise direction and 129 points in the 
cross-section. This structure permits the removal of sections along the wing where other topologies can be 
used to enhance grid resolution. For each change in topology a decoupling of grid density is introduced. 
Since a wing is typically modeled with a C-grid topology, the middle section of the O-grid that extends 
along the side of the shuttle is replaced with a C-grid. Additionally the tail grid is spanwise mapped with a 
C-grid. The C-grid on the wing has an additional 48 cells on the upper and lower surfaces for a total of 129 
points on the wing, and 225 points in the cross-section that includes the fuselage. The C-grid on the tail 
places 49 points vertically along the leading edge, without adding more points in the streamwise direction. 
Thus, enhanced grid resolution of the wing surfaces is enabled through the addition of different topologies. 

When the damage scenarios are modeled, this region is replaced with an O-grid because grid density 
from the edges of the C-grid that connect to the edges of the damage geometry is independent of the 
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(a) O-grid. 


(b) C-grid. 


(c) H-grid. 


Figure 2: Basic structured grid topologies on a wing. 


remaining sections of the volume grid. 

The remaining dimension in this volume grid is the body to outer boundary direction, which was 
modeled with 65 points. However, if further resolution is required in the boundary layer or bow shock 
region, this dimension could be easily increased because it is consistent among all the blocks of the 
common baseline volume grid. 

Thus, using different topologies to model different geometric features of the Shuttle Orbiter while 
maintaining a consistent computational direction for the boundary layer, grid density was enhanced 
without requiring additional points in regions were the grid resolution was already adequate for the 
accurate capture of flow field phenomenon by the CFD codes. 

GRID QUALITY 


Most flow simulation software will work with any structured grid, providing the computed cell volumes are 
always positive. In order to obtain “accurate” flow simulations, however, very specific conditions of the 
grid must be met. As noted in the section above, there are physical characteristics which must exist, such 
as orthogonality. The most common method of determining the applicability of a grid towards accurate 
flow simulations is to measure the physical characteristics, also known as grid metrics. 

The grid quality metrics used to assess the grids generated for use by the ACAI team were based on 
grid generation experience at NASA Langley Research Center (LaRC). These grid metrics include 
cell-to-cell stretching, interior and near wall orthogonality, and cell volume. Based on previous work in 
structured grid generation, grid stretching should be less than a factor of 1.5[5, 6]. During the CAI 
analyses a more stringent requirement of monotonic grid stretching variation was required by the USA 
software to obtain accurate heating analyses. Grid line intersection within the volume and near the wall of 
the vehicle should be less than 10 degrees from an orthogonal condition[7], for the grid to be considered to 
be of high quality. 

These metrics were used to assess the grid quality by evaluating the maximum stretching, minimum 
orthogonality, and root mean square (RMS) values of the orthogonality and stretching. The RMS values 
were evaluated in combination with a factor of three times the standard deviation to identify the limits of 
the metrics. Most grids that did not adhere to these metrics were either re-constructed or modified until 
the metric limits were satisfied. Based on anecdotal evidence obtained during the use of the grids 
generated for this investigation, failure to adhere to limits of these metrics can result in inaccurate flow 
simulations and could slow the convergence of flow field residuals with the CFD software. 
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THE PROCESS 


The build-up of structured grid topologies shown in figure 1 was chosen for the following reasons: 

1. There would be no singularity in the flow domain. 

2. The embedded O-grid would permit the local re-construction of a section of the grid to accommodate 
related damage scenarios, without requiring full configuration regeneration. 

3. The embedded O-grid topology focuses grid resolution to the vicinity of the damage where high 
grid-point resolution is required. 

4. Grid adaptation in the OML surface normal direction remains an option for those solvers that can 
use it. 

5. Piecewise construction techniques could be used to spread the workload among different participants. 

The first and third issues address the accuracy of the flow simulation software, as previously identified. 

The topology choice was made to permit rapid changing of the shuttle orbiter geometry, leading towards a 
flexible construction process and flexible usage of geometry for complex flow field simulation and analyses. 
Flexibility in geometric usage is derived from the ability to simply swap volume grid components in order 
to build different shuttle orbiter configurations. Each of the damage scenarios that were examined by the 
CAI team required the meshing of an interior cavity region in addition to the exterior mesh previously 
discussed. The interior of the cavity region is fitted with either an H-grid or a double O-grid, depending on 
the damage scenario and flow solver that was in use. 

The construction techniques that are enabled through the choice of these exterior and interior 
topologies are as varied and different as the participants on the ACAI team. Within the team, JSC 
provided direction for the types of geometries to be modeled, while ARC and LaRC shared the 
responsibilities of constructing the volume grids, utilizing the process identified in figure 3. In figure 3 the 
steps that ARC were responsible for performing are highlighted in blue, while the steps that LaRC 
performed are shown in red, and combined responsibilities of ARC and LaRC are shown in purple. Each 
volume grid to be generated used the same process identified in figure 3. The wall surface grid that was 
used to model the geometry of the Shuttle Orbiter Columbia was generated by ARC. This usually is one of 
the most time consuming steps in the grid generation process because the grid is restricted to the OML of 
the orbiter. All grid generation and smoothing of the wall grid is performed on the computer aided design 
(CAD) definition of the vehicle. Conversely, the interfaces of blocks, and the other faces that define a 3D 
block are not as restricted. Thus, generating 3D surfaces in general 3D space as well as the volume grid 
that fills each block requires approximately the same amount of work as constructing the surface grid. 

During the use of the grid generation process, several hurdles had to be overcome, including: 

1. Obtaining an outer mold line (OML) definition of the shuttle orbiter. 

2. Generating block interfaces that are orthogonal to the wall of the vehicle. 

3. Maintaining C-2 (curvature) continuity across matching block boundaries. 

4. Limiting grid stretching and skewness at the surface and within the volume. 

Each of these issues are addressed in the following sections. Contained in these sections are also 
explanations of how the grid generation process unfolded and how the delegation of responsibilities lead to 
the use of LaRC and ARC resources. This process lead to the successful completion of structured grids for 
viscous flow simulations in a rapidly developed and high quality product. 


ORBITER OUTER MOLD LINES 

At the time of the Columbia accident, the most accurate CAD definition of the vehicle OML that was 
available was that developed by Mississippi State University under contract to JSC in the mid 1990’s. 
Inspection of the surface geometry identified several flaws in wing shape and elevon trailing edge. Hence, a 
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Figure 3: Grid generation process. 
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new OML was sought. JSC, who had already been working on an improved vehicle OML, provided this 
new orbiter definition within the first few weeks of the Columbia accident. 

Due to some geometry nuances and minor flaws in the new OML from JSC, an effort was undertaken 
by both ARC and LaRC to construct a “watertight” OML of the shuttle orbiter using two different 
methods. The method utilized by ARC was designed to obtain an OML shape that could be used to assess 
aerothermodynamic issues. Specifically, ARC appended small regions of false surfaces to the Orbital 
Maneuvering System (OMS) pods and the wing to simplify the geometry, as shown in figure 4. These 
appendages are outside of the domain of influence. Breaklines along the actual geometry were maintained 
so that subsequent aerodynamic force integration could be performed on solution data obtained during the 
aerothermodynamic analyses. Addition of the geometry, though small, provided the necessary OML for 
structured surface and volume grid generation. 



Figure 4: Ames Research Center OML modifications. 

The second watertight OML developed by LaRC, resulted in a more accurate representation of the 
orbiter needed for aerodynamic simulations. This OML surface was not used for the ACAI team because 
the details were too difficult to accommodate with a structured grid. No superfluous surfaces were added, 
and strong surface discontinuities present in the actual Shuttle OML were maintained. Each OML was 
used in the investigation, and differed only on the top and trailing edge of the vehicle. Since 
aerothermodynamic characteristics are not usually influenced by these regions, the differences in geometry 
were considered negligible. 


GENERATING HIGH QUALITY GRIDS 


With the OML established for aerothermodynamic simulations, the process identified in figure 3 
commenced. As noted in figure 1, the topology of the grid was agreed upon prior to generating any grids. 
During the investigation it became clear that several grids would have to be generated. Although the 
topology provides a working structure of where damage scenarios to the orbiter could be easily 
accommodated, the detailed local regriding required to support each specific damage scenario was still 
necessary. The GridGen[8] software was used to generate surface grids. Several versions of the wall surface 
grid were evaluated by LaRC for skewness and stretching during the wall surface grid generation process. 
The Tecplot[9] software in conjunction with the CFD Analyzer[10] (CFDA) was used to perform the 
computations for evaluating grid skewness and stretching. The CFDA was chosen because of its integration 
into the Tecplot environment and because it was able to assess grid stretching across matching block 
boundaries. Each consecutive grid generated had consistently increasing quality, as the skewness and 
stretching were minimized. 
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BLOCK INTERFACE ORTHOGONALITY ENHANCEMENT 


Once the wall surface grid was complete, the volume grid construction commenced. This portion of the 
process was performed by ARC and LaRC in concert with one another. Due to the construction techniques 
available, controlling the block interfaces at the wall of the orbiter was very difficult. Although the 
GridGen software does have the ability to project grid lines orthogonal to the wall, generating the entire 
interface with an edge along the wall surface such that the interior grid of the interface is orthogonal to the 
wall, is impossible. There exists no method within GridGen to ensure that the interface grid along the wall 
of the Orbiter is orthogonal. Hence, an alternative technique had to be used. 

Through extensive work at LaRC on the X-37 program, a method to improve an existing surface grid 
for orthogonality was developed. This technique, as applied to the CAI, requires an initial surface grid and 
the vectors to control the location of the grid along an edge. The vectors, typically obtained from the 
CFDA, are combined with the initial surface grid in the Volume Grid Manipulator [11] (VGM) tool to 
produce an orthogonally corrected interface. Briefly, referring to figure 5 for vector construction and 
figure 6 for the corrected grid, the method entails the following steps: 

1. Input the initial surface grid and the controlling vectors along an edge. 

2. Reverse the controlling vectors along an edge. 

3. Construct a second set of vectors interior to the interface grid, but close to the edge where corrections 
are to be performed. 

4. Utilize Hermite Vector Interpolation[12] (HVI) to regenerate a region close to the edge in question by 
interpolating the vector fields. 

5. Perform Trans-Finite Interpolation[13] on a small region near the ends of the modified edge grid. 

6. Output the corrected grid for elliptic smoothing. 



Figure 5: Vectors used within VGM to correct orthogonality of a block interface grid near the wall. 

The resulting grid, as shown in figure 6 is orthogonal to the wall based on the adherence to the vector 
field used to create it. This VGM created surface grid is then used as a defining surface within GridGen, to 
construct a new surface mesh. A final surface grid is created by elliptically smoothing the new surface 
mesh within GridGen using both wall orthogonality and a fixed surface condition to ensure both 
orthogonality to the OML and orthogonality to the surface grid edges. Each interface connected to the wall 
is corrected in the same manner with the VGM software. This correction ensures that the volume grid 
based on these interfaces can be orthogonal to the wall of the orbiter, thus satisfying a desired grid 
attribute for the CFD solvers. Additionally, by manipulating the grid in VGM, the process is fast and 
efficient, so rapid generation of the volume grid is not impeded. 




Figure 6: Corrected block interface grid near the wall using HVI and TFI. 

IMPOSING C-2 CONTINUITY WHILE 
LIMITING SKEWNESS AND STRETCHING 


Upon completion of generating and correcting the interfaces of the volume grid to be constructed, the 
volume grid is generated. Volume grid generation for the Shuttle Orbiter was performed primarily using 
GridGen, and 3DGRAPE/AL[14] where GridGen failed. In particular, the topology of the tail region was 
implemented using complex surface grids which made it extremely difficult for GridGen to generate a grid 
without crossed grid lines. As such, this region was generated with the 3DGRAPE/AL software. The most 
difficult issue to overcome in the volume grid generation was the construction of continuous grid lines 
across matching block boundaries. Typically, software such as GridGen uses an orthogonal grid condition 
across such a boundary. Although this does produce a continuous grid line across the matching block 
interface, the resulting grid is difficult to control on the interior of the domain, and inflections in the grid 
can lead to excessive skewness of the mesh within the volume domain. There are several methods for 
producing slope continuous grid lines across matching block boundaries, but each has disadvantages that 
are not conducive to rapid structured volume grid generation. For example, use of a Neumann boundary 
condition or averaging the point locations at the first cell away from the interface will change the interface 
entirely. Another technique has been to swap the vectors from one block into the other that share an 
interface. The latter method is used by GridGen and can require multiple exchanges of angles across a 
matching block interface before the slope continuous grid line is produced. Regardless of the method used, 
to date there has been no boundary condition that will preserve the interface and provide slope continuous 
grid lines across matching block interfaces for rapid volume grid generation. 

To overcome the problems associated with the development of slope continuous grid lines across 
matching block boundaries, a new technique to correct the grid was developed specifically for the CAI. The 
method makes use of vector field averaging. The steps involved include: 

1. Computing vector fields from the first cell away from the interface on both sides or in each block that 
share the interface. 

2. Reversing one vector field from one side of the interface. 

3. Averaging the vectors at the interface. 

4. Reversing a copy of the vectors to the side where the vectors from step 2 originated. 

5. Using the vectors to project the first cell points to correct the first layer of grid away from the 
interface in both blocks that share the interface. 

The steps above perform the correction. However, to use the correction effectively, the correction has 
to be performed several times. As shown in figure 7a with a representative plane from a volume grid on the 
Orbiter, the grid is initially generated with 3DTFI. Then the correction is applied as shown in figure 7b, 
and the volume grid is smoothed again, but maintaining the angles and cell sizes produced from the 
correction at the interface, resulting in the grid shown in figure 7c. This cycle of smoothing and correcting 
is performed 2 to 3 times in order to obtain the slope continuous condition across the matching block 
interface. The result of using this method is a high fidelity structured volume grid with slope continuous 
grid lines across all matching block interfaces. 
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Figure 7: A single cycle for applying the matching block interface continuity correction. 

VOLUME GRIDS PRODUCED 

Utilization of the grid generation process as identified above with the new technology and the teamwork, 
resulted in a rapid generation of volume grids. The first mesh was generated using the flexible framework 
process for the baseline undamaged configuration shown in figure 1. The vast majority of problems related 
to mesh defects such as grid stretching and grid line intersection skewness, and data transfer problems 
between LaRC and ARC were discussed and resolved using this first mesh. For example, the initial OML 
surface mesh generated by ARC and analyzed by LaRC had poor stretching and boundary matching 
parameters. Subsequent iterations between ARC and LaRC resolved these issues. 

Upon completion of the baseline common grid, the process had matured to the extent that subsequent 
grids where damage was included, were generated efficiently with very few iterations required to obtain the 
desirable quality. To evaluate flow conditions and phenomenon as related to the CAI, four full body 
exterior damage scenario volume grids were generated using the established process in figure 3. Each grid 
was generated by removing the yellow region in figure 1, and regenerating the region to reflect the different 
damage scenarios, including those shown in figure 8. 

In each of the geometries shown in figure 8, cavities were generated and inserted where the red 
topology lines are drawn. For the RCC Panel 6 and 9 cases, multiple cavity grid topologies were developed 
for half and full panel missing configurations, but each used the same exterior surface and volume mesh. In 
each cavity case the inner face of the exterior volume mesh in the region directly above the intended cavity, 
was developed to correspond to the original OML surface. Thus, each of these damage scenario grids could 
be used without an additional internal cavity grid, as a surrogate for calculating a smooth, undamaged, 
OML case. The cavity grids were mated to these exterior meshes along the cavity outer boundary/ OML 
surface. In most cases the interior cavity grids were constructed by the individual organizations responsible 
for calculating the particular damage scenario case. All cavities, except RCC-8, were modeled with various 
boundary conditions to simulate a pressurization of the wing leading edge. Only RCC-8 had an interior 
volume grid that extended beyond the cavity itself and included a portion of the channel of the wing 
leading edge. By limiting the scope of the volume grid generation to the specific damage scenario location, 
resources were focused on the generation of high fidelity grids for the region of the wing leading edge. This 
significantly reduced the time required to generate a new full body volume grid because the wing leading 
edge region was simply swapped for the different grids. Note that in some cases, the grid density was 
modified along the fuselage blocks to enable finer grid resolution in the cavities to be modeled. 

Nonetheless, the process used and the limiting of scope of the grid generation enabled rapid volume grid 
generation for the CAL 


SUMMARY 


A flexible framework for the generation of high fidelity structured volume grids for use with the Columbia 
Accident Investigation (CAI) was developed. The framework consisted of a grid generation process and a 
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(a) Main landing gear (MLG). 



(c) Reinforced carbon-carbon (RCC) panel 8. 

Figure 8: Four damage scenario surface geometries 
of the CAI. 


(b) Reinforced carbon-carbon (RCC) panel 6. 



(d) Reinforced carbon-carbon (RCC) panel 9. 
topologies for the aerothermodynamic simulations 


division of labor among the aerothermodynamic CAI team members. The process was a compromise given 
the the aggressive time line and wide range of flow solver capabilities. Given more time, other best 
practices in grid generation could have been utilized. The topology utilized in this framework for damage 
scenario grids was an embedded O-grid, which enables grid metric customization with respect to the 
geometry alteration without adversely affecting grid requirements on the remainder of the vehicle. The 
embedded O-grid topology also permitted the construction of multiple geometric changes through the 
regeneration of a limited region of the baseline volume grid. To further reduce the grid density 
requirements on the grid, a C-grid topology was used on the wing and tail components. The grid generation 
process utilized the strengths of decoupled surface and volume grid generation along with evaluation of grid 
quality to ensure that the final grid for computational use would be acceptable to all participants in the 
aerothermodynamic computations of the CAI. Utilizing a combination of teamwork between NASA Centers 
and the grid generation process resulted in a framework for rapid generation of high fidelity grids. 
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